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Abstract: This paper presents an approach for realizing a phase difference measurement of 
a new gyro. A silicon micromachined gyro was mounted on rotating aircraft for aircraft 
attitude control. Aircraft spin drives the silicon pendulum of a gyro rotating at a high speed 
so that it can sense the transverse angular velocity of the rotating aircraft based on the 
gyroscopic precession principle when the aircraft has transverse rotation. In applications of 
the rotating aircraft single channel control system, such as damping in the attitude 
stabilization loop, the gyro signal must be kept in sync with the control signal. Therefore, 
the phase difference between both signals needs to be measured accurately. Considering 
that phase difference is mainly produced by both the micromachined part and the signal 
conditioning circuit, a mathematical model has been established and analyzed to determine 
the gyro's phase frequency characteristics. On the basis of theoretical analysis, a dynamic 
simulation has been done for a case where the spin frequency is 15 Hz. Experimental 
results with the proposed measurement method applied to a silicon micromachined gyro 
driven by a rotating aircraft demonstrate that it is effective in practical applications. 
Measured curve and numerical analysis of phase frequency characteristic are in accordance, 
and the error between measurement and simulation is only 5.3%. 
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1. Introduction 

With the increasing development of MEMS and inertial guidance technology, all kinds of 
micromachined gyros have successfully developed and are gaining increasing popularity for shared use 
in military and civil applications [1-3]. According to the driving structure, a MEMS gyro can be 
divided into two types. One is a gyro with a driving structure, and another is gyros without driving 
structures. The vast majority of reported micromachined rate gyroscopes utilize a vibratory proof mass 
suspended by flexible beams above a substrate. The primary working principle is to form a vibratory 
drive oscillator, coupled to an orthogonal sense accelerometer by the Coriolis force. Using anchor, the 
proof mass is suspended above the substrate, making the mass free to oscillate in two orthogonal 
directions-the drive and the sense directions (Figure 1) [4-8]. 



This kind of gyro is difficult to design and manufacture and the cost is high. In order to avoid the 
difficulties brought about by the driving part, this paper puts forward a novel gyro that uses the 
rotation of the aircraft itself as a driving part. In this paper, the studied gyro, which has been fabricated 
on s monocrystalline silicon wafer by means of bulk micromachining manufacturability techniques, 
belongs to the gyro without driving structure class. Since there is no driving structure, the structure is 
simple and easy to process [9,10]. This new gyro is used on rotating aircraft for flight attitude control. 
The gyro sensing element is called a silicon pendulum. Aircraft spin provides angular momentum for 
the silicon pendulum so that it can sense the transverse angular velocity when aircraft appears to turn. 
Using detection circuits, the angular vibration of the silicon pendulum is transformed into an 
alternating electric signal. The signal frequency is same as the spin frequency of aircraft, and the signal 
envelope is proportional to the transverse angular velocity. 

In the practical applications, dynamic characteristics of micromechanical gyroscope is the key [11-17]. 
For a single channel control of aircraft, the novel gyro is for damping in the whole stable loop, as 
shown in Figure 2. In Figure 2, Uk is the control signal and Ug is the gyro's output signal. The term Us is 
the control signal of the steering engine and it equals the value of Uk minus Ug. The term Us controls the 



Figure 1. A generic MEMS gyroscope with the driving part. 
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aircraft rotation via the actuating mechanism and the gyro senses a transverse angular velocity of f2 , 
then the gyro outputs a signal Ug and feeds it back to the end of the control signal, so a stable control 
loop can be obtained based on the novel gyro in the single channel control system of the rotating 
aircraft. The measurement object of the gyro is the input transverse angular velocity of the rotating 
aircraft, where Ug is the gyro's output signal such that Ug = KQcos((p - fi), where K is scale factor, O 
and g) are respectively the input angular velocity and spin angular velocity of the aircraft, and [5 is the 
phase difference between Uk and Ug. 

Figure 2. The block diagram of damping loop. 
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Because the gyro signal is a negative feedback, the gyro signal and control signal need to be in 
phase. Therefore, it is a key that the phase difference be accurately measured so as to adjust and 
compensate the phase in practical applications. 

2. Kinematic Equation 

2.7. Kinematic Equation 

Figure 3 shows the gyro principle structure diagram, where the middle section is a sensitive element, 
called a silicon pendulum. Through bulk micro-mechanical processing technology, it is fabricated on a 
mono-Si wafer. There are two copper plated ceramic electrodea on both sides of the silicon pendulum, 
on which copper has been plated by mean of a sputtering process. The plate electrode and silicon 
pendulum form symmetrical detection capacitances, where m is a detection capacitance positive pole 
and n is another detection capacitance positive pole, as shown in Figure 4. 



Figure 3. Gyro structure schematic diagram. 



elastic torsional beam 



copper plated ceramic 
electrode 




silicon pendulum 



Sensors 2013, 13 



11054 



Figure 4. Detection capacitances. 
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Figure 3 shows that elastic torsion beams are located on both ends of the constraint center of the 
silicon pendulum and they are fixed on the frame. Nitrogen is encapsulated between the electrode and 
the silicon pendulum. When the silicon pendulum is vibrating around the constraint center, the 
capacitance changes periodically. Then, by picking up the capacitance circuit, the gyro outputs a 
sensitive signal. 

The rotating aircraft-driven silicon micromachined gyro is installed on the rotating aircraft. In 
Figure 3, 6>y is the output axis, called also precession axis, ox is the input axis and oz is the driving 
axis, which is accordance with the spin axes of a rotating aircraft. When it is rotating at a high speed 
around the spin axis, the aircraft will drive the silicon pendulum rotating at the same angular velocity 
as the aircraft spin angular velocity (p and the silicon pendulum will also acquire angular momentum. 

Angular momentum is same in the direction as the oz axis. If the aircraft turns with angular velocity 
of Q around the input axis, the silicon pendulum will be forced to produce precession. Thus, along 
with the direction of precession oy axis, the precession moment acts upon silicon pendulum. The 
precession moment can be balanced with various moments such as the inertia moment, damping torque 
and elastic moment. Concerning these moments, because the silicon pendulum is vibrating around the 
precession axis, the angular acceleration generates inertia moment. Having angular velocity, the 
squeeze film resistance of nitrogen generates a damping torque, and the elastic torsion beam also 
generates elastic moment. In short, when aircraft undertakes a forced turn around the input axis ox, as a 
result, the silicon pendulum will generate precession to make angular momentum tend to the input 
transverse angular velocity Q. The term d represents the angular velocity of the silicon pendulum 
vibration. Suppose that the input angular velocity i2 is a constant, with the aircraft spin, the silicon 
pendulum will generate a simple harmonic vibration. After the detection circuit and conditioning 
circuit, the gyro will output a periodic voltage signal, as shown in Figure 5. A case is introduced to 
validate the correctness and applicability of the Si micromachined gyro driven by a rotating aircraft. 
The experimental platform is shown in Figure 6. 

Figure 6 shows that the gyro was installed on an emulator and emulator was mounted on an angular 
vibration table. The spin angular velocity was given by the emulator and the value of the spin angular 
velocity had been adjusted to 5,4007s (15 Hz). Input transverse angular velocity was provided by the 
angular vibration table and was inputted a constant value of 180°/s. The gyro's output signal was 
measured via a scope and the measured waveform is shown in Figure 7. If the aircraft is not rotating at 
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a constant rate but rather a sinusoidal function with a frequency of 1 Hz, the gyro's output signal is a 
AM signal and the envelope of the signal corresponds to the sinusoidal function, as shown in Figure 8. 
The test control panel is shown in Figure 9. 

Figure 5. Output voltage waveform with a constant input angular velocity. 
U 




Figure 6. Experimental platform. 
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Figure 8. AM waveform. 
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Figure 9. Test control panel. 




Waveforms in Figures 5 and 7 are perfectly in accordance. Therefore, testing validates that the 
previous analysis is correct. With respecting to the above-mentioned motion, it can be illustrated by 
rigid motion around fixed point. Therefore, using the Euler dynamic equation based on rigid 
movement around a fixed point and its projection on the procession axis, we will easily derive the 
motion equation to express Equation (1) as in [10]: 

Jyd + Dd + \^J^ - J + Kj^a = {^J J ^^f2(pcos{^(pt) (1) 

where J ^ , J ^ and J ^ represent the moment of inertia of the silicon pendulum relative to x, y and z, 
respectively. Kj^ is the torsion stiffness of the elastic torsion beam and £) is a damping factor, d 
represents the value of the vibration angular velocity of the silicon pendulum, f2 represents the value 
of the input angular velocity, and (p represents the value of the spin angular velocity of the aircraft. To 
simplify Equation (1) further, it can be written as: 

d+2gco^d + co^a = /q cos {cpt) (2) 



where: 
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(3) 



Equation (2) is a second order constant coefficient linear differential equation, then, the solution of 
Equation (2) is given by: 

a=A^e~"' cos^^cOyt + S^ + A^cos^gX-/^) (4) 

where Aj is a integration constant, which value is determined by the initial conditions, n is an 
attenuation factor, including n = go)^ and oj^ = '>^o)^'^-r^ . Equation (4) indicates that the solution 
consists of a damping vibration and a harmonic vibration. The first item of Equation (4) would be 
quickly attenuated to zero, then the forced vibration will reach steady state and it can be calculated 
as follows: 



a = 



^[(j^-j^-j^y+K,J+{D^y 

Then, the amphtude of the angular vibration can be written as: 

(j^+j^-j,)<p 



tCOS 



(5) 



(6) 



The phase difference is: 



P = arctan 



^ 2n<p ^ 



(7) 



2.2. Kinetic Parameters 



Figure 10 is the structural diagram of the silicon pendulum. The gyro frame is fixed during 
packaging and cannot move. There is a movable part, called the silicon pendulum, inside the frame. Its 
two ends are elastic torsion beams, which connect with the frame. There are 14 damping bars on two 
sides of the silicon pendulum. Squeeze film damping, has a big impact on the dynamic response of 
MEMS microstructure, and can be used to adjust the quality factor of the micro-mechanical structure. 
Increased damping bars can reduce the pressure of the gas film, and the gas film damping is reduced 
accordingly, so in the MEMS design, in order to speed up the release or decrease damping, we usually 
include a release bar or so-called damping bar in the silicon pendulum. 
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Figure 10. Silicon pendulum structure. 
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In the chip process, we have four steps. The first step shapes a movable section, the second step 
etches the damping bars, the next is buffer layer etching on the beams and the last step forms the 
silicon pendulum. The masks of the four steps and chip are shown in Figure 11. 



Figurell. The masks of four steps. 




Using the above four step etching, the final chip is obtained. A picture of the chip is shown in 
Figure 12. 

Figure 12. Picture of the chip. 
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2.2.1. Torsion Stiffness 

Considering that it is convenient to calculate the torsion stiffness of an elastic torsion beam, it is 
supposed that on processing the torsion angle is proportional to the length of the beam, warping of the 
cross sections of elastic torsion beam are same, and values are same, but in opposite direction, for the 
torsion moment of the two ends of an elastic torsion beam. 

Under the condition of the above assumptions, using elastic mechanics, the torsion stiffness with 
rectangular cross section can be obtained by: 



512G^'w ^ 1 



^ ^ «=1,3,5,... n 



2t , riTrw 

1 tann 

riTTW 2t 



(8) 



where L , w , and t are the length, width and thickness of the elastic torsion beam, respectively. G is 
the sheer modulus of the silicon material. 

Practically the torsion angle a is limited to 0.003 rad, is large enough, thus we can simplify 
Equation (8), and then Equation (9) can be obtained: 

^ 2 Gt'w 

If G , L , w , and t are chosen to be: 

L = 600 |Lim, w = 400 |Lim, f = 75 |Lim, G = 5. IxlO'' nW (10) 

the calculated result is: 

TT^ =2.507x10"' N-m (11) 



2.2.2. Damping Factor 



The displacement of the silicon pendulum edge node is far less than both the lateral size and gap 
relative to the electrodes in the torsion motion of Figure 4. By using numerical FEA simulation, the 
maximum displacement is only 886 nm. Therefore, in the analysis of the problem, we make the 
assumption that the silicon pendulum undergoes a rigid motion of small amplitude, the gas film gap is 
just a time function, and the linear solution of the Reynolds equation is given for calculating the 
damping coefficient of the silicon pendulum. 

When a rectangular plate with length A and width B is moving relative to the undersurface of a gap 
distance of h , the damping coefficient of the squeeze film resistance can be expressed by: 



f = 



damp 

dh/dt 



AB'ju 



1- 



192 

Att' 



^ — tank 



n=l,3,5,.. 



~2B 



(12) 



where // is the gas viscosity coefficient. 

Since the infinite series can quickly converge, the damping coefficient is approximately equal to the 
first term of Equation (12), that is: 



/= 



1 192 , 7v' 

1 tann — 

Att^ IB 



(13) 
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In order to calculate the damping factor D shown in Figure 10, the first step is to figure out the 
damping factor of the entire silicon pendulum, then, to calculate the damping factor of the damping 
hole, notch and fourteen damping bars. Finally, the entire silicon pendulum damping factor minus the 
damping factors of the missing parts above. As a result, the real damping factor of the silicon 
pendulum is obtained. 

The chosen relevant dimensions are: h = 20 |im, A = B =14 mm, ju = 17.81 x 10~^ Pa-s, thus, the 
damping factor of the silicon pendulum can be figured out: 

D = 1.676xlO~^N-m-s (14) 



2.2.3. Moment of Inertia 

Moment of inertia ( , , ) relative to the principal axes of inertia are calculated in the same 
way as the damping factor. Finally, the obtained moment of inertia can respectively be written as: 



= 2.27794 g-mm^ 
73; = 2.1527 g-mm^ 
7^ = 4.43064 g-mm^ 



(15) 



3. Phase-Frequency Characteristics 

3.1. Phase-Frequency Characteristics 



In order to analyze the phase-frequency characteristics of the Si micromachined gyro driven by a 
rotating aircraft, the vibration of the silicon pendulum can be represented as a SISO system, of which 
the excitation signal is given by x[t) and the response signal is given by j(^). The term x[t) will 
simulate the excitation generated by the moment of procession. The term y(^t) simulates the 
displacement of the angular vibration of the silicon pendulum. Thus, Equation (2) can be expressed in 
another form as follows: 



Using the Laplace transfonnation, the transfer function of Equation (16) can be written as: 

1 



(16) 



H{s) = 



(17) 



Then, the frequency characteristic is: 



H{j(o) = - 



^ 1 ^ 



1- 



r \2 
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f \ 

CO 



(18) 
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Finally, we can easily obtain the phase-frequency characteristic: 



2? 



CO 



cp{^a)) = arctan 




(19) 



In actual applications, the working frequency (co) is less than the natural frequency of co^ , so 
Equation (19) shows that the output signal has a phase lag relative to the input, or said differently, the 
silicon pendulum angular vibration and input angular velocity have a phase delay. When the spin 
angular velocity of the aircraft reaches the natural frequency of the silicon pendulum, this phase lag is 
90 degrees. 

3.2. Phase-Frequency Characteristics of Dijferent Damping Ratios 

The phase-frequency characteristics are simulated and plotted for a Si micromachined gyro driven 
by a rotating aircraft. The plotted curve is shown in Figure 13. 



In Figure 13, it is clearly shown that when the damping ratio is smaller, the phase difference also 
becomes less. Practically we can always design to an appropriate damping ratio, such as 0.5 or 0.7. 

3.3. Conditioning Circuits Impact 

The Si micromachined gyro principle block diagram is shown in Figure 14, which shows that the 
phase difference is not only relevant to the vibration of the silicon pendulum, but also affected by the 
conditioning circuits. Phase difference is equal to the sum of the phase difference of two parts. The 
sconditioning circuits consist of a differential amplifier circuit, band-pass filter circuit, low-pass filter 
circuit and gain circuit. 



Figure 13. Phase frequency characteristic curves for different damping ratios. 
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Figure 14. Block diagram of gyro working principle. 
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We chose an AD620 chip, which constitutes the differential amplifier circuit. We connect the gyro's 
sensed signal with its noninverting and inverting inputs. Its magnification times are adjusted by a 
resistor cross input. Meanwhile, we add a LPF at the noninverting and inverting inputs, respectively. 
The transfer function of the differential amplifier circuit is written as: 

25.7 



10"'^ + 1 



(20) 



Because the gyro's working frequency ranges from 1 Hz to 50 Hz, the band width of the BPF is 
designed to be 81 Hz. The transfer function for BPF is: 

0.64^^+3520^ 



H,{s) = 



(21) 



0.64/ +324^ + 2000 

In order to filter the noise, we design the second order LPF. Its cut-off frequency is 63.5 Hz and 
transfer function is: 

1 



H,{s) = 



6.618x10"'/ +3.52xlO"'^ + l 



,-3 . 



(22) 



For the gain circuit, its function is mainly amplifying the signal and large transmission bands. The 
peak gain is 17.3 dB, so the transfer function is designed to be: 

118.252^ + 68000 



HAs) = — 

' 4.96x10^/ +16.213S + 22000 

Then, for the entire conditioning circuits, the transfer function is: 

H{s) = H,{s)H^{s)H,{s)H,{s) 



(23) 



(24) 



According to Equation (24), we solve the amphtude-frequency characteristic and phase frequency 
characteristic. The Bode plot is as shown in Figure 15. 

According to the amplitude-frequency characteristics of Figure 15, in the range of l'-50 Hz, the 
gain is relatively flat. In addition, the phase frequency characteristics of Figure 15 show that the phase 
difference curve is approximately linear. 
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Figure 15. Bode plot of Equation (24). 
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4. Dynamic Simulation of Phase Difference 

4.1 . Simulation Modeling 

The model structure built according to the working principles of the gyro is shown in Figure 16. It 
includes the second order vibration model of a silicon pendulum and the conditioning circuit model. 



Figure 16. Modeling of the dynamic simulation. 
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In Figure 16, an input periodic signal simulates the excitation of the second-order vibration system. 
The frequency of the periodic signal is the spin frequency of the aircraft. The constant transverse 
angular velocity is modulated on the sinusoidal wave amplitude. The second order vibration system is 
a closed loop system. Due to the vibration angular velocity of the silicon pendulum, the damping 
moment, which is formed by the damping factor D , feeds back to the input. Meanwhile, the torsion 
moment, which is formed by the torsion stiffness Kj. , also feeds back to the input. To summarize, the 
precession moment, damping moment and torsion moment together form a loop with negative 
feedback. There are two differential operators S in the loop model. The conditioning circuit models 
are located on the last step of the modeling. 

4.2. Dynamic Simulation 

The driving signal frequency is set to 15 Hz and simulates the aircraft spin frequency. Because the 
input drive signal simulates the aircraft rotation driving, when the aircraft has a transverse angular 
velocity, the silicon pendulum will generate an angular vibration. Because the silicon pendulum is of 
microscopic scale, the moment of inertia around the driving axis is very small, so the angle of vibration 
is also very small and the detected signal is a small signal. Therefore, the input driving signal 
amplitude is smaller in the modeling. In order to compare the phase difference between the input 
driving signal and the gyro output signal, a gain of between the driving signal and the oscilloscope 
is added for phase difference comparison. Using dynamic simulation, the simulation results are shown 
in Figure 17. 

Figure 17. Dynamic simulation waveform. 
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Figure 17 indicates that the gyro output frequency is same as the driving signal and their 
frequencies are 15 Hz. The phase delay is 8 ms, that is, the phase lag is: 

^(p^lnfM (25) 

n 

where / is the spin frequency of the aircraft, is the phase delay. If F = 75 Hz, z/^ = 5 ms are 
chosen then the calculated phase difference isAg) = 43.2°. 
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5. Measuring of Phase Difference 

5.1. Measuring Platform 

The gyro was installed on a MEMS three-axis precision turntable. Figure 18 shows the installation 
on the turntable. Figure 19 is the control table, and Figure 20 shows the working status. 

Figure 18. Measuring platform. 
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Figure 19. Control platform. 




Figure 20. Working status. 
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Figure 20 shows that the roUing axis was rotating at a certain speed and generated angular 
momentum for the sihcon pendulum. When the input transverse angular velocity was generated by the 
yaw axis or pitch axis, the gyro senses an input transverse angular velocity. In order to measure 
the phase difference with respect to different frequencies, we adjust different rolling velocities via the 
control platform in Figure 19. Then we compare the phase between the gyro's output signal and 
angular displacement of the rolling axis, and we obtain the phase difference with respect to different 
frequencies. 

5.2. Measurements 



Given an input angular velocity of 180°/s, the turntable rolling axis is rotating at a high speed. At 
the same time, through the data acquisition card of the MEMS three-axis precision turntable, we will 
obtain data of the frame angular displacement and gyro output signal. After data acquisition is 
completed, the speed of the rolling axis is increased. In the same way, measurement is accomplished. 

In order to compare to the simulation result as shown in Figure 17, we choose sampled data 
of (p = 15 Hz, Q = 7507s. Using these sampled data, we plot a curve, as shown in Figure 21. The blue 
line represents the turntable rolling axis angular displacement, while the red one represents the gyro 
output signal. 

Figure 21. Phase delay comparison chart. 




As shown in the Figure 21, the phase difference between the gyro output signal and driving reaches 
a value of 40.9°. Under the different rolling frequency, the calculated phase differences are recorded in 
Table 1. The fitting curve of the measured phase frequency characteristic is shown in Figure 22. 



Table 1. Phase-frequency characteristic experimental records. 
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Figure 22. The experimental phase-frequency characteristic curve. 
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6. Results and Discussion 

Figure 22 shows that the characteristic phase frequency curve, ranging from 5 Hz to 25 Hz, is in 
accordance with the simulation curve in Figure 13. Therefore, the results demonstrate that the 
measurement method is correct, but because of the maximum limits of the rolling frequency of the 
turntable rolling axis, the range of measured frequency is not outside the limit of 30 Hz. For dynamic 
measurements of a spin frequency of 15 Hz, the simulation result of the phase difference is 43.2° and 
the measured value is 40.9°. The error is only 5.3%, so the measured phase difference can meet the 
needs of adjusting the phase in single channel control applications of rotating aircraft. 

7. Conclusions 

This paper introduces a new type of micromachined gyro, and its working principles were analyzed. 
For the MEMS micromachined gyro driven by a rotating aircraft, the movement equation is presented 
and the dynamic parameters are calculated. Based on the motion equation, we analyzed the phase 
frequency characteristics and plotted the characteristic phase frequency curve under different damping 
ratio scenarios. Combined with the dynamic parameters, a dynamic simulation model has been set up. 
When the spin frequency is 15 Hz, the obtained phase difference simulation is 43.2°. On the basis of 
this analysis, we designes a phase delay measurement method. Using this method, we have measured 
the characteristic phase frequency curve. Ranging from 5 Hz to 25 Hz, the measured curve was in 
accord with the simulation. When the spin frequency is chosen as 15 Hz, by measurements we 
obtained a phase difference of 40.9°. The relative error with respect to the simulation is only 5.3%. 
The measured results prove that the accuracy is good and the data of the measured phase difference 
can meet the requirements of adjusting the phase for the single channel control of a rotating aircraft. 
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